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We show a geometric formulation for qubit state discrimination, that can be generally applied, 
and provide the complete solution in a closed form when arbitrary qubit states are given with equal 

a priori probabihtics. It is shown that the guessing probability does not depend on detailed relations 
among given states, such as angles between them, but on a property that can be assigned by the set 
of given states itself. This also shows how a set of quantum states can be modified such that they 
give the same guessing probability. The general structure of optimal measurements is characterized, 
which also explains that no measurement sometimes gives an optimal strategy. 
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Discrimination of quantum states is a fundamental 
processing to extract information encrypted in collected 
quantum states. In practical applications, its framework 
charac;tc;rizes communication capabilities of encoding and 
decoding messages via quantum states [1] [2]. Its useful- 
ness as a theoretical tool to investigate quantum informa- 
tion theory has also been shown in secure communication 
[3] or randomness extraction [4]. 

Despite of much effort so far [5], however, apart from 
two-state cases [1], the discrimination task is generally 
put into a numerical procedure [6] and approximated. 
Otherwise, restricted cases when some specific symmetry 
is given can be analyzed, e.g. Ref. [7]. In fact, little 
is known in general about optimal discrimination once 
an arbitrary set of quantum states is given. Needless to 
speak about the importance in its own right, the lack of a 
general method for state discrimination even in simple in- 
stances is, due to its fundamental importance, obviously 
and also potentially a significant obstacle preventing fur- 
ther investigation in both quantum information theory 
and quantum foundation. 

The present work makes remarkable development in 
the problem of state discrimination, in particular, for 
cases of qubit states, and leads to significant improve- 
ment along the line. We analyze optimality conditions 
of the minimum-error discrimination, and provide a ge- 
ometric formulation of optimal discrimination of qubit 
states. We show the guessing probability, i.e. the max- 
imal probability of making correct guesses, in a closed 
form when they are given with equal a priori probabil- 
ities, and present a method of finding optimal discrimi- 
nation using the geometry of qubit states. Optimal mea- 
surements are characterized accordingly. The discrimi- 
nation is shown for arbitrary N qubit states in which 
no symmetry may exist. From these, it is shown that 
the guessing probability does not depend on detailed re- 



lations among given states in general but the property 
assigned by the set of given states. This shows how a 
set of states can be modified such that the modification 
cannot be recognized in the discrimination task. 

For the purpose, let us briefly summarize the 
minimum-error discrimination in the context of a com- 
munication scenario of two parties, Alice and Bob. They 
have agreed with N alphabets {x}^;^ and states {pxj^Li, 
as well as a priori probabilities {qx}x=i- Alice's encod- 
ing works by mapping alphabet x to state px, and re- 
lating states with a priori probabilities {Qx}^i- This 
can be seen as Alice's pressing button x with probability 
gfx, and then Bob's guessing among states {pxIxLi given 
with a priori probabilities {^xj^i, which we write by 

{9x,Px}^=l- 

Bob's discrimination of quantum states is described by 
Positive-Operator- Valued-Measure (POVM) elements: 
{Mx > 0}^^i satisfying J^^M^ = I. Let PB|A(x|y) de- 
note the probability that Bob has a detection event on 
Mx that leads to the conclusion that Px is given, while a 
state Py is provided by Alice's sending message y. This 
is computed as follows, PB|^(x|y) := P(x|y) = tr[Ma;Py]. 
The figure of merit is the maximal probability that Bob 
makes a correct guess on average, called the guessing 
probability, 

-Pguess = niax y'qxtr[Mx/9x], y'Mx = /, (1) 

L xJx— 1 X X 

where the maximization runs over all POVM elements. 
This naturally introduces the discrimination as an opti- 
mization task. 

In fact, the above can be put into the framework of 
semidefinite programming [6]. A useful property in this 
approach is that a given maximization (minimization) 
problem can be alternatively described by its dual, a min- 
imization (maximization) problem. The dual problem to 
the maximization in Eq. (1) is obtained as follows. 
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Pguess = mintr[i4:], K > q^p^, Vx = 1, • • • ,N. (2) 
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In this case, the minimization works to find a single pa- 
rameter K which then gives the guessing probabihty, e.g. 
the approach in Ref. [8]. 

For convenience, wc call the problem in Eq. (1) as 
the primal, with respect to the dual in Eq. (2). Note 
that solutions of two problems do not coincide in gen- 
eral. The fact that the guessing probability can be ob- 
tained from both primal and dual optimizations follows 
from the property called strong duality. This holds when 
both primal and dual problems have a non-empty set of 
parameters satisfying given constraints, which is referred 
to as the feasible problems. Once both problems are fea- 
sible, the strong duality holds, and then it follows that 
solutions from both problems coincide each other. 

Apart from solving those optimization problems, there 
is another approach called a complem,entarity problem. 
This collects optimality conditions that parameters of 
both primal and dual problems have to satisfy to give 
optimal solutions. Then, any set of parameters satisfy- 
ing optimality conditions immediately provide optimal 
solutions of primal and dual problems. As more param- 
eters arc taken into account, this is not considered to be 
easier, however, the advantage lies at its usefulness to 
find general structures of a given optimization problem. 

In the semidcfinitc programming formulation, the op- 
timality conditions can be summarized by the so-called 
Karus-Khun- Tucker (KKT) conditions. For quantum 
state discrimination, they are given by, together with two 
constraints in Eqs. (1) and (2), 

K = q^p^ + rx(Tj;, and (3) 
Mr[axMJ =0, Vx = I,.- - ,iV (4) 

for a set of complementary states {(Jx}x=i with non- 
negative coefficients {r^ > 0}x=i POVM elements 
{Mx}^i. Once states {rx,<Tx}^i and measurements 
satisfying these conditions are found, they arc automat- 
ically optimal to give solutions in both primal and dual 
problems. Prom the fact that the strong duality holds 
in this case, it is clear that the guessing probability is 
obtained from either problem. Note that the former con- 
dition in Eq. (3) is called the Lagrangian stability, and 
shows the existence of a single operator K that can be 
decomposed into N different ways. The latter one in 
Eq. (4) is the complementary slackness that shows the 
orthogonality relation between primal and dual parame- 
ters. 

The particular usefulness of KKT conditions in state 

discrimination is that, as it is shown in the above, they 
separate the guessing probability (i.e. tr[K] see Eq. (2)) 
from optimal measurements: a single operator K solely 
characterizes the guessing probability and optimal mea- 
surements itself are independently expressed. The opera- 
tor K can be explained to have A'' decompositions of q-x/Ox 
and TxCTx for each x. Optimal mcasiircments can be de- 
scribed as POVM elements orthogonal to states {<Jx}^=i 



for each x. The discrimination problem is then equivalent 
to finding states {crxj^i to fulfill these conditions. 

We now show a geometric formulation to find comple- 
mentary states. Let us first define the polytopc of given 
states {(?xiPx}^i denoted as 7^({(?xi Px}^i) in the un- 
derlying state space, in which each vertex corresponds to 
Qxpx- It is useful to rewrite the condition in Eq. (3) as, 

fePx - QyPy = TyUy - TxPx, V X, y (5) 

This shows that two polytopes 'P({<Zx, Px}^i) of given 
states and P({rx, <Tx}^i) of complementary states to 
search for arc actually congruent. Thus, the striicture of 
complementary states is already determined from given 
states {<7x,/5x}^i- Once the state geometry is clear e.g. 
[9], the formulation can be applied. 

For qubit states, their geometry can generally be de- 
scribed in the Bloch sphere using the distance measure, 
Hilbert- Schmidt norm. In what follows, wc restrict our 
consideration to qubit states and apply the geometric 
formulation to discrimination among them. For a qubit 
state /?x, we write its Bloch vector as v(/9x), so that 
Px = (-^ + v(px) • (?)/2 where <? = (X, F, Z) Pauli ma- 
trices X, F, and Z. 

For qubit states, we can characterize the general form 
of optimal measurements from the KKT condition in Eq. 
(4). Suppose that > 0, otherwise, the measurement 
can be arbitrarily chosen. To fulfill the condition, it is not 
difficult to sec that optimal POVM elements are either of 
rank-one [10] or the null operator. If CTx = |V'x)(V'x| then 
Mx = m^ltp^) {tp^\ with coefficients m^, where it holds 
that v(^x) = — v('!/'x )• If ^ state cTx is not of rank-one, the 
only possibility to fulfill the KKT condition in Eq. (4) is 
that the measurement corresponds to the null operator, 
i.e. Mx = 0. In fact, no measurement can sometimes 
give optimal discrimination [11]. Note that, however, 
as measurements arc done in most cases (otherwise, the 
completeness of POVM elements in the following is not 
fulfilled), one does not have to immediately assume that 
{<^x}^=i are not of rank-one from the beginning. Then, 
for cases where measurements are done, corresponding 
complementary states must be of rank-one - otherwise, 
the orthogonality in Eq. (4) cannot be fulfilled. 

Once states {crx}x=i are found, projectors for opti- 
mal measurement are determined accordingly. What re- 
mains is that POVM elements fulfill the completeness, 
-^x = I, or equivalently in terms of Bloch vectors, 
X^x'^xv(V'x) = while X^2,mx = 2 with rrix > 0, Va;. 
This refers to finding a convex combination {m^}^^-^ of 
Bloch vectors {v(f/'x such that it results to zero, 

i.e. the origin of the Bloch sphere. This is equivalent 
to the condition that the convex hull of Bloch vectors 
{v('(/'x)}^i of complementary states contain the origin 
of Bloch sphere. As we will show later, this is always 
fulfilled by complementary states. To summarize, once 
complementary states are found, it is automatic to have 
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optimal POVMs as their Bloch vectors are determined 
and the completeness is also straightforward. 

We can thus proceed to construction of comple- 
mentary states in the Bloch sphere for qubit states 
{9x,Px}x=i- Let us identify the polytope 7'({gx, Px}^=i) 
in the state space as the convex hull of their Bloch 
vectors {(jr^, v(px)}^i, so that each vertex ^xPx corre- 
sponds to the Bloch vector (/xv(px)- Then, the task is 
to find the polytope ^({rx, a-x}^i) of complementary 
states that is congruent to ViiqxT Px}x=i) ^^"^ Bloch 
sphere. Moreover, as it is shown, most of complementary 
states {(Tx}^i are of rank-one and themselves lie at the 
border, thus, the polytope P{{<t^}^^i) only of them is 
maximal in the Bloch sphere. Using these, a geometric 
approach can be generally employed. 

All these already give the guessing probability in a sim- 
ple way for cases when qubit states are given with equal 
a priori probabilities, that is, when ^x = 1/-^ for all x. 
In this case, it is not difficult to see a general form of 
the guessing probability. Substituting Qx = 1/-/V in the 
KKT condition in Eq. (3), it is obtained that parameters 
{rx}^i are equal. We put r := for all x. This holds 
true for arbitrary set of quantum states in general. Then, 
the guessing probability is written as 

Pguess = t,[K] = 1 + r, with r = (6) 

iV ||(Tx — CTy II 

where the equation for the parameter r is from the rela- 
tion in Eq. (5), and the distance measure can be taken 
as the Hilbert-Schmidt norm Dus (as it is natural in the 
Bloch sphere) or the trace norm Dt- Both measures give 
the same value of r since they are related only by a con- 
stant: Diis = V^Dt for qubit states. This means that, 
the parameter r can be obtained by either of distance 
measures while referring to the geometry in the Bloch 
sphere, since the parameter is only a ratio, see Eq. (6). 

The parameter r corresponds to the ratio between two 
polytopes, V{{l/N, px}^^i) of given states and the other 
■PdcTx}^!) only of complementary states, as it is shown 
in Eq. (6). We recall that most of {(t}^i are pure (i.e. 
rank-one), lying at the border of the Bloch sphere. This 
implies that the polytope 'P({o'x}^i) of complementary 
states is clearly the maximal in the Bloch sphere. In this 
way, the polytope 'P{{ax}x=i) always contain the ori- 
gin of the sphere, from which optimal measurements can 
be constructed. Finally, two polytopes Pdl/TV, pxl^i) 
and 'P{{<Tx}^^i) are similar from the relation in Eq. (5) 
by the ratio r, since two polytopes '^({1/.^, Px}^i) and 
■p({l/A^, (Txl^i) are congruent. 

We summarize the method of discrimination. 

1: Construct a polytope from given states as the convex 
hull of {1/N,v{px)}x=i ill Bloch sphere where 
vertices are \{p-x)/N. 

2: Expand the polytope such that it keeps simiar to the 
original one and is also maximal in the Bloch sphere 




FIG. 1. (A) Six states {1/6, /9x}x=i in the half-plane are given 
with purities {/x}x=i, i.e.. OPx ~ /x/6, where three of them 
(x = 1, 3, 4) have the same purity and the others are less. See 
also that the relation in Eq. (5) is fulfilled. The ratio r in Eq. 
(6) can be obtained by expanding the given polytope until it 
is maximal in the plane. This is also the ratio between radii 
of two circles covering respective polygons. Thus, Pgucss = 
l/6-f/i/6. Complementary state CTx corresponds to OSx. For 
X = 2, 5, 6, they are not pure states and thus no measurement 
on these states is optimal [10]. (B) For {l/A'', pxj^Li of pure 
states, each vertex of the polyhedron corresponds to the Bloch 
vector of state px/N. The ratio r is equal to Q1Q2/R1R2, see 
also Eq. (5), and thus Pgucss = 2/N. Even if these states 
are modified, if the minimal sphere covering the polyhedron 
is unchanged, the guessing probability remains the same. 

(as most of {cx}^i are pure), and then compute 
the ratio r of the resulting polytope with respect to 
the original one. The guessing probability is thus 

obtained, Pgucss = ^/N + r. 

3: Rotate the maximal polytope within the Bloch sphere 
until it is found that corresponding lines are paral- 
lel to the original ones to fulfill Eq. (5). From this, 
corresponding vertices are complementary states 
W^}x=i and optimal POVM elements are explic- 
itly constructed accordingly. 

This completely solves the problem of discrimination of 
qubit states given with equal a priori probabilities. 

It is already observed that the guessing probability 
does not depend on detailed relations of given quantum 
states to discriminate among, but a property from the 
whole set {px}^i, since its ratio r is the relevant param- 
eter. If given states are modified such that the polytope 
has the same ratio r, then the same guessing probability 
is given. This means that, in the communication sce- 
nario we have introduced in the beginning, Alice who 
encodes messages can choose, or modify, sets of quan- 
tum states {pxj^Li in such a way that Bob who decodes 
from quantum states cannot recognize her modification 
using the optimal guessing. This actually defines equiva- 
lence classes of sets of quantum states in terms of optimal 
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guessing [12]. 

In the following, we apply the method to various cases 
of qubit state discrimination. The simplest example, also 
the case when a general solution is known, is for N = 2, 
say pi and p2- Following the instruction in the above, i) 
the polytope constructed by given two states corresponds 
to a line connecting two Bloch vectors of the states. The 
length can be computed using the trace distance as — 
P2II/2. Then, ii) the maximal ploytope similar to the 
original one is clearly the diameter of the Bloch sphere, 
which has length 2 in terms of the trace distance, hence, 
r = \\pi — P2II/4 (which equivalently can be obtained 
with the Hilbert-Schmidt distance). Substituting this in 
Eq. (6), the Helstrom bound in Ref, [1] is reproduced. 
Then, iii) the diameter can be rotated until it is parallel 
to the original one. Thus, optimal measurements are also 
obtained. 

Next, let us consider N states on the half- 
plane. We can begin with those states {1/N, px}^^i 
which are equally distributed and generally not pure. 
They are characterized by Bloch vectors: y{Px) = 
/x(cos ^x, sin 6'x, 0) where 6^ = 27rx/iV. For these, no 
general solution is known except cases called geometri- 
cally uniform states, that {/x}^i are equal [7]. Here, 
suppose that N is even so that we can also assume 
/jv/2 = fN = niaxx/x- Then, applying the method in- 
troduced, one can easily find that the ratio depends on 
the maximal purity, that is, r = f^/N, and the guessing 
probability is obtained as Pguess = ^/N + Jn /N , no mat- 
ter what purities the other A''— 2 states have, see also Fig. 
1. The result in Ref. [7] is also reproduced. The assump- 
tion of the equal distribution on angles can be relaxed 
while keeping On =-n + 9^/2 and fM/2 = In = maxx /x, 
for which the guessing probability is also the same no 
matter how other N — 2 states are structured. This is 
because, as it is shown in Eq. (6), they are given with 
equal probabilities and the ratio r is unchanged. 

Optimal measurements can be analyzed as follows, 
based on the geometric formulation, see also Fig. 1. 
For two states Pn/2 and pj^ having the maximal pu- 
rity, it is clear that measurement is applied, and let 
o-jv = \ilJN){i'N\ and (7jv/2 = |'^i'iv/2)(V'JV/2|- Prom these, 
optimal POVM elements are straightforward. For the 
other states, say {p^} having fz < Jn, it holds from the 
KKT condition in Eq. (5) that vjv — Vz is parallel to 
r(v(crAr) — v(crz)). This simply shows that their comple- 
mentary states {(Tz} cannot be pure states, i.e. not of 
rank-one. Then, corresponding POVM element is the 
null operator, that is, no measurement on these states 
gives an optimal strategy. 

The method can be applied to a set of qubit states 
having a volume. For instance, let us look at the case 
when pure states are given such that their Bloch vec- 
tors form a regular polyhedron of N vertices, see Fig. 
1. Following the instruction in the above, the parame- 
ter r can be obtained as the ratio of two spheres, one 



the Bloch sphere and the other the minimal sphere cov- 
ering the polyhedron of given states {1/N, pxj^Li- From 
this, we have r = 1/N, and the guessing probability is 
thus, Pgucss = 2/N. One can also modify angles between 
those N states such that the minimal sphere covering 
the polyhedron remains the same, and then the guessing 
probability is unchanged 

The presented method can be in principle applied to, 
high dimensional states once their geometry is clear, or 
qubit states with unequal a priori probabilities. For the 
latter, although the geometry is clear, we do not have yet 
a general and systematic method to derive the guessing 
probability. For these cases, interesting and illuminating 
examples are to be presented elsewhere [12]. 

To conclude, we have shown a geometric formulation 
for qubit state discrimination and provide the guessing 
probability in a closed form for equal a priori proba- 
bilities. This makes a significant contribution to the 
study of quantum state discrimination. Optimal mea- 
surements are characterized accordingly. It is shown how 
qubit states can be modified while the guessing proba- 
bility remains the same. As qubits are units of quantum 
information processing, we envisage that the method of 
discrimination and results presented here would be useful 
to develop further invc;stigations of qubit applications, or 
approaches to related open questions e.g. Ref. [13]. 
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